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Abstract 

' y-H ' ■ It is known that entangled mixed states that are positive under partial transposition (PPT states) 

must have rank at least four. In a previous paper we presented a classification of rank four entan- 
gled PPT states which we believe to be complete. In the present paper we continue our investiga- 
tions of the low rank entangled PPT states. We use perturbation theory in order to construct rank 
five entangled PPT states close to the known rank four states, and in order to compute dimensions 
and study the geometry of surfaces of low rank PPT states. We exploit the close connection be- 
tween low rank PPT states and product vectors. In particular, we show how to reconstruct a PPT 
state from a sufficient number of product vectors in its kernel. It may seem surprising that the 
number of product vectors needed may be smaller than the dimension of the kernel. 
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in ; 1 Introduction 

T— I ! 

t^j- \ Quantum entanglement between subsystems of a composite physical system is a phenomenon which 

clearly distinguishes quantum physics from classical physics Q. Entangled quantum states show 
correlations between measurements on the subsystems which can not be modelled within classical 
physics with local interactions. A classical model would have to be a joint probability distribution 
of quantitities that are incompatible in the quantum theory, and the existence of such a joint proba- 
bility distribution, consistent with locality, implies so called Bell inequalities Q, or even equalities 
as in the three particle states known as GHZ states, introduced by Mermin, Greenberger, Home, and 
Zeilinger @1. The correlations in entangled quantum states violate Bell inequalities and GHZ 
equalities. 

A pure classical state of a composite system has no correlations between measurements on subsys- 
tems, since classical measurements are deterministic. A statistical ensemble of pure classical states, 
what we may call a mixed classical state, can have correlations, but these correlations can not violate 
Bell inequalities, by definition. 

The only pure quantum states that are not entangled are the pure product states, which resemble 
pure classical states in that they have no correlations at all. By definition, a mixed quantum state is a 
statistical ensemble of pure quantum states, and it is said to be separable if it can be mixed entirely 
from pure product states. The separable mixed states are not entangled, since they can not violate Bell 
inequalities. The entangled mixed states are precisely those that are non-separable. For this reason, 



1 



the mathematical distinction between separable and non-separable mixed states is important from the 
physical point of view. 

The separability problem, how to characterize the set S of separable mixed states and decide 
whether a given mixed state is separable or not, is known to be a difficult mathematical problem Q. 
It motivates our work presented here and in previous papers, although we have not studied so much 
the separable states directly as the larger class of mixed states called PPT states |6l[7j[8j[9]. 

The separable mixed states have the property that they remain positive after partial transposition, 
they are PPT states, for short. The set V of PPT states is in general larger than the set S of separable 
states, but the difference between the two sets is surprisingly small in low dimensions, and in the very 
lowest dimensions, 2 x 2, 2 x 3, and 3x2, there is no difference ifTDI . 

The condition of positive partial transpose is known as the Peres separability criterion liTTI . It is a 
powerful separability test, especially in low dimensions where the difference between the two sets V 
and S is small. It can be used for example to prove that any pure quantum state is either entangled or 
a pure product state. 

We study especially the lowest rank entangled PPT states, on the assumption that they are the 
easiest ones to understand. In the present paper we discuss in particular how to construct PPT states 
of rank 4 and 5 in 3 x 3 dimensions, and rank 6 in 4 x 4 dimensions. A central theme is how the states 
are constrained by the existence of product vectors in the kernel. Another central theme is perturbation 
theory, which we use to conduct rank 5 PPT states close to rank 4 states, and to study surfaces of PPT 
states of fixed low rank. We compute numerically the dimensions of such surfaces, and we show how 
to follow a surface by numerical integration of an equation of motion. 

The relation between PPT states and product vectors 

The close connection between PPT states and product vectors has been used earlier, for example to 
prove the separability of sufficiently low rank PPT states lfl2Tl . 

Bennett et al. lfT3l[T4l introduced a method for constructing low rank mixed states that are obvi- 
ously entangled PPT states, using what they called Unextendible Product Bases (UPBs). From a UPB, 
defined as a maximal set of orthogonal product vectors which is not a complete basis of the Hilbert 
space, one constructs an orthogonal projection Q and the complementary projection P = t — Q. Then 
p = P/(Tr P) is an entangled PPT state. 

The UPB construction is most successful in the special case of rank 4 PPT states in 3 x 3 di- 
mensions. In Ref. 13 we argued, partly based on evidence from numerical studies, that an extended 
version of the UPB construction, including nonunitary but nonsingular product transformations on the 
states, is general enough to produce all rank 4 entangled PPT states in 3 x 3 dimensions. 

Unfortunately, attempts to apply the UPB method directly in higher dimensions fail, even when 
the kernel contains product vectors, because there can not exist a sufficient number of orthogonal 
product vectors. The orthogonality is essential in the construction by Bennett et al. of the PPT state as 
a projection operator. We would like to generalize the construction in such a way that it works without 
the orthogonality condition. 

One possible generalization is to construct projection operators as more general convex combina- 
tions, or even as non-convex linear combinations, of pure product states. This idea is explored in a 
separate paper [15 1. 

In the present paper we discuss in general the constraints imposed on a PPT state p by the existence 
of product vectors in its kernel, and we show that these constraints are so strong that they actually 
determine the state uniquely. A surprising discovery is that in cases where the kernel contains a finite 
overcomplete set of product vectors, the state p can be reconstructed from only a subset of the product 
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vectors, and the number of product vectors needed may even be smaller than the dimension of the 
kernel. 

From this point of view, the important question is exactly what conditions the product vectors 
must satisfy in order for the constraint equations to have a solution for p. We can answer this question 
in the familiar special case of rank 4 PPT states in 3 x 3 dimensions, but not in other cases. We 
consider this an interesting problem for future research. 

Outline of the paper 

The contents of the present paper are organized as follows. First we review some linear algebra, in 
particular degenerate perturbation theory, in Sections El and |4] The main purpose is to introduce 
notation and collect formulas for later reference. 

In Section[5]we discuss the rank 4 PPT states in 3 x 3 dimensions. We review the UPB construction, 
based on orthogonal product vectors in the kernel, before we describe an approach which is different 
in that the product vectors need not be orthogonal. The new approach also throws some new light on 
a set of reality conditions that limit the selection of product vectors to be used for constructing rank 4 
PPT states. 

In Section [6] we discuss the rank 5 PPT states in 3 x 3 dimensions. We find an 8 dimensional 
surface of rank 5 PPT states in every generic 5 dimensional subspace, but we have not found any 
general method to construct such states. However, we show how to construct rank 5 PPT states by 
perturbing rank 4 PPT states. Again, the product vectors in the kernel of the rank 4 state play an 
important role in our construction of the rank 5 states. 

In Section [7] we discuss rank 6 PPT states in 4 x 4 dimensions. The kernel of such a state has 
dimension 10, and contains 20 product vectors. The remarkable result we find is that the state can be 
constructed from only 7 product vectors in the kernel. An arbitrary set of 7 product vectors does not 
produce a rank 6 PPT state, but we do not know how to select sets of product vectors that can be used 
in such a construction. 

In Section [8] we discuss briefly how to determine numerically the dimensions of surfaces of PPT 
states of fixed rank. We find that the dimensions are given by a simple counting of independent 
constraints, except for the very lowest rank states, for which the constraints are not independent. 

Finally, we discuss in Section [9] how to study a surface of PPT states by numerical integration of 
equations of motion for curves on the surface. In this way one may study for example the curvature 
of the surface, or how a curve on the surface approaches the boundary of the surface. 

2 Some basic linear algebra 

2.1 Density matrices 

Let Hn be the set of Hermitean N x N matrices. It has a natural structure as a real Hilbert space of 
dimension N 2 with the scalar product 

(X,Y) = Tt(XY) . (1) 
A mixed state, or density matrix, is a positive Hermitean matrix of unit trace. We define 

V = V N = {peH N \p>0,Trp = l} . (2) 
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Because it is Hermitean, a density matrix p has a spectral representation in terms of a complete set of 
orthonormal eigenvectors ipi € <C N with real eigenvalues Aj, 

N N 

P = ^2 Xi With = S V ' 1 = ' ^ 

i=l i=l 

The rank of p is the number of eigenvalues Aj 7^ 0. The pseudoinverse of p is defined as 

P + = . w 

it is equal to the inverse p~ x if p is invertible. The matrices 

P = p + P = PP + = ]T ^ , Q = 1 - P = titi (5) 

i.Aj^O i,Aj=0 

are Hermitean and project orthogonally onto two complementary orthogonal subspaces of (D , P 
onto Img p, the image of p, and Q onto Ker p, the kernel of p. The relations Pp = pP = PpP = p 
and Qp = pQ = QpQ = will be used in the following. 

We say that p is positive, or positive semidefinite, and we write p > 0, when Aj > for i = 
1,2, . . . , N. An equivalent condition is that ip*pifi > for all ift G C . It follows from the last 
inequality and the spectral representation of p that tp^px/j = if and only if pip = 0. 

The definition of positive Hermitean matrices by inequalities of the form ip* pip > implies that 
T> is a convex set. That is, if p is a convex combination of p\,p2 € V, 

p = Api + (1 - A)p 2 with 0< A < 1 , (6) 

then p € V. Furthermore, since Ker p = {tp \ tp^ pip = 0} when p > 0, it follows that 

Ker p = Ker p\ n Ker p2 , (7) 

independent of A, when p is a convex combination as above. Since Ker p is independent of A, so is 
Imgp = (Kerp)- 1 . 

A convex set is defined by its extremal points: those points that are not convex combinations of 
other points. The extremal points of V are the pure states of the form p = ipip^ with ip G C . Thus, 
the spectral representation in eq. d3) is an expansion of p as a convex combination of N or fewer 
extremal points of V. 



Finite perturbations 

In the following, let p be a density matrix and define the projections P and Q = 1 — P as in eq. (f5]). 
Consider a perturbation 

p' = P + eA, (8) 

where A ^ is Hermitean, and Tr A = so that Tr p' = Tr p. The real parameter e may be finite or 
infinitesimal, we will first consider the case when e is finite. 

We observe that if Img A C Img p, or equivalently if PAP = A, then there will be a finite range 
of values of e, say ei < e < €2 with ei < < €2, such that p' € V and Imgp' = Imgp. This 
is so because the eigenvectors of p with zero eigenvalue will remain eigenvectors of p' with zero 
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eigenvalue, and all the positive eigenvalues of p will change continuously with e into eigenvalues of 
p>. 

The other way around, if p' € V for e\ < e < €2 with e\ < < £2, then p' is a convex combination 
of p + e\A and p + €2^4 for every e in the open interval e\ < e < 62- Hence Img p' is independent of 
e in this interval, implying that Img A C Img p and PAP = A. 

This shows that p is extremal in V if and only if there exists no A ^ with TV A = and 
P^4P = A. Another formulation of the condition is that there exists no p' € V with p' ^ p and 
Img p' = Img p. A third equivalent formulation of the extremality condition is that the equation 

PAP = A (9) 

for the Hermitean matrix A has A = p as its only solution (up to proportionality). In fact, if PBP = B 
and Tr B ^ 0, then we have PAP = A and Tr A = when we take 

A = B-(TvB)p. (10) 

Infinitesimal perturbations 

Assume now that Img A (£. Img p. The question is how an infinitesimal perturbation affects the zero 
eigenvalues of p. When p is of low rank we need degenerate perturbation theory, which is well known 
from any textbook on quantum mechanics. 

To first order in e, the zero eigenvalues of p are perturbed into eigenvalues of p' that are e times 
the eigenvalues of QAQ on the subspace Kerp. Similarly, to first order in e, the positive eigenvalues 
of p are perturbed into positive eigenvalues of p', in a way which is determined by how p and PAP 
act on Img p. 

It is clear from this that, to first order in e, the condition 

QAQ = (11) 

is necessary and sufficient to ensure that the rank of p' equals the rank of p, and that p' > both for 
e > and for e < 0. 

More generally, to first order in e, the rank of p' equals the rank of p plus the rank of QAQ. For 
example, if we want to perturb p in such a way that the rank increases by one, then we have to choose 
A such that 

QAQ = a #t , (12) 

where <j) £ Ker p is a normalized eigenvector of QAQ with a ^ as eigenvalue. Since QAQ is 
Hermitean, a must be real. If a > 0, then p' > for e > but not for e < 0. 

Projection operators on 

Using the projections P and Q defined above we define projection operators on Hn, the real Hilbert 
space of Hermitean N x N matrices, as follows, 

PI = PXP , 

QX = QXQ = X - PX - XP + PXP , (13) 
RX = (I - P - Q)X = PX + XP - 2PXP . 
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Here I is the identity on H^. It is straightforward to verify that these are complementary projections, 
with P 2 = P, Q 2 = Q, PQ = QP = 0, and so on. They are symmetric with respect to the natural 
scalar product on Hn, for example, 

(X, PY) = Tr(XPYP) = Tr(PXPY) = (PX, Y) . (14) 

Hence they project orthogonally, and relative to an orthonormal basis of they are represented by 
symmetric matrices. 

Relative to an orthonormal basis of <D N with the first basis vectors in Img p and the last basis 
vectors in Ker p, a Hermitean matrix X takes the block form 

w)> (15) 

with W = U and = W. In this basis we have 

«-(!?). 

and hence, 



2.2 Composite systems 
Product vectors 

If N = NaNb then the tensor product spaces = <S) (a complex tensor product) and 
Hn = Hn a <8> Hn b (a real tensor product) may describe a composite quantum system with two 
subsystems A and B of Hilbert space dimensions Na and Nb- 
A vector ip € C w then has components ipj = tpij, where 

7=1,2,. ..,N o ij = 11, 12,..., 1AT B , 21, 22,..., iV A iV B . (18) 

A product vector = 4> ® \ has components ^ij = 0iXj- We see that -0 is a product vector if and 
only if its components satisfy the quadratic equations 

Y^-Yto - 4>uiPkj = o . (19) 

These equations are not all independent, the number of independent complex equations is 

K = (N A -1)(N B -1) = N-N A -N B + 1. (20) 

For example, if V'n / Owe get a complete set of independent equations by taking i = j = 1 and 
k = 2, 3, . . . , N A , I = 2, 3, . . . , N B . 

Since the equations are homogeneous, any solution tp / gives rise to a one parameter family of 
solutions ctp with c G C A vector ip in a subspace of dimension n has n independent complex com- 
ponents. Since the most general nonzero solution must contain at least one free complex parameter, 
we conclude that a generic subspace of dimension n will contain nonzero product vectors if and only 
if 

n>K + l = N-N A -N B + 2. (21) 
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The limiting dimension 



n = N - N A - N B + 2 (22) 

is particularly interesting. In this special case a nonzero solution will contain exactly one free param- 
eter, which has to be a complex normalization constant. 

Thus, up to proportionality there will exist a finite set of product vectors in a generic subspace of 
dimension n = N — N A — N B + 2. The number of product vectors is |H] 

= ( N A + N B -2\ (N A + N B -2)\ 
P V *A-1 J (N A -1)\(N B -1)\ ■ { } 

A generic subspace of lower dimension will contain no nonzero product vector, whereas any subspace 
of higher dimension will contain a continuous infinity of different product vectors (different in the 
sense that they are not proportional). 



Partial transposition 

The following relation between matrix elements, 

defines X p , the partial transpose of the Hermitean matrix X with respect to the second subsystem. 
A density matrix p is called separable if it is a convex combination of tensor product matrices, 

9 = ^2 P k ° k ® Tk > ( 25 ) 

k 

with cjfc € T>n a , Tfc £ T>N B ,Pk > 0, YlkPk = !• We denote by S the set of separable matrices. 
The partial transpose of the above separable matrix is 

p p = Y, pk(7k ®( T ^ T - ° ■ (26) 

k 

The positivity of p p is known as the Peres criterion, it is an easily testable necessary condition for sep- 
arability. For this reason it is of interest to study the set of PPT (Positive Partial Transpose) matrices, 
defined as 

V = { P \ P p >o} = T>nv p . (27) 

We may call it the Peres set. A well known result is that V = S for N = N A N B < 6, whereas V is 
strictly larger than S in higher dimensions ifTOl . 

We will classify low rank PPT states by the ranks (m, n) of p and p p , respectively. Note that ranks 
(m, n) and (n, m) are equivalent for the purpose of classification, because of the symmetric roles of 
p and p p . 
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Product transformations 

A product transformation of the form 

p^t fj = aVpV ] with V = V A ® V B , (28) 

where a is a normalization factor and V A € SL(A^4,C), Vb € SL(iVs,C), preserves positivity, 
rank, separability, and other interesting properties that the density matrix p may have. For example, it 
preserves positivity of the partial transpose, because 

(p' f = aVp p V ] with V = V A ® Vg . (29) 

The image and kernel of p and p p transform in the following ways, 

Img// = V Imgp , Kerp' = (V^)' 1 Kerp , (30) 

and 

Img(p') P = V\mgp p , Ker(p') P = (^ f ) _1 Kerp p . (31) 

All these transformations are of product form and hence preserve the number of product vectors in a 
subspace. 

We say that two density matrices p and p' related in this way are SL <S> SL equivalent, or simply 
SL equivalent. The concept of SL equivalence is important to us here because it simplifies very much 
our efforts to classify the low rank PPT states. 

3 Restricted perturbations 

We have seen that eq. © ensures that the perturbation p' = p + eA preserves the image of p, so that 
Img p' = Img p for infinitesimal values of the perturbation parameter e, and Img p' C Img p even for 
finite values of e. The weaker condition in eq. dTTT> ensures only that the rank of p' equals the rank of 
p for infinitesimal values of e. 

We want to discuss how to use perturbations with similar restrictions in order to study, for example, 
the extremal points of the convex set V. In particular, we are interested in perturbations that either 
preserve the ranks (m, n) of p, or else change these ranks in controlled ways. 

In a similar way as we did for p, we define P and Q = 1 — P as the orthogonal projections onto 
Img p p and Ker p p . Then we define 

FX = (PX P P) P , 

QX = (QX P Q) P = X- (PX P ) P - {X P P) P + (PX P P) P , (32) 
Rl = (I — P — Q)X = {PX P ) P + (X P P) P - 2{PX P P) P . 

These are again projections on the real Hilbert space Hn, like P, Q and R, again symmetric with 
respect to the natural scalar product on H^. 

We may now use the projection operators on Hjy to impose various restrictions on the perturbation 
matrix A. 
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Testing for extremality in V 

The extremality condition for V is derived in a similar way as the extremality condition for V based 
on eq. ©. Clearly p is extremal in V if and only if there exists no p' € V, p' ^ p, with both 
Img p' = Img p and Img(p') p = Img p p . Another way to formulate this condition is that A = p is 
the only solution of the two equations PA = A and PA = A. 

Since P and P are projections, the equations PA = A and PA = A together are equivalent to 
the single eigenvalue equation 

(P + P)A = 2A . (33) 

They are also equivalent to either one of the eigenvalue equations 

PPP.4 = A , PPP A = A . (34) 

Note that the operators P + P, PPP, and PPP are all real symmetric and therefore have complete 
sets of real eigenvalues and eigenvectors. In fact, the eigenvalues are all non-negative, because the 
operators are positive semidefinite. 

When we diagonalize P + P we will always find A = p as an eigenvector with eigenvalue 2. If it 
is the only solution of eq. (l33l . this proves that p is extremal in V. If A is a solution not proportional to 
p, then we may impose the condition Tr A = (replace A by A — (Tr A)p if necessary), and we know 
that there exists a finite range of both positive and negative values of e such that p' = p + eA £ V, 
hence p is not extremal. 

It should be noted that in our numerical calculations we may find eigenvalues of P + P that differ 
from 2 by less than one per cent. However, since eigenvalues are calculated with a precision close 
to the internal precision of the computer, which is of order 1CP 16 , there is never any ambiguity as to 
whether an eigenvalue is equal to 2 or strictly smaller than 2. 

Perturbations preserving the PPT property and ranks 

The rank and positivity of p is preserved by the perturbation, to first order in e, both for e > and 
e < 0, if and only if QA = 0. Similarly, the rank and positivity of p p is preserved if and only if 
QA = 0. These two equations together are equivalent to the single eigenvalue equation 

(Q + Q)A = . (35) 

Again Q + Q is real symmetric and has a complete set of real eigenvalues and eigenvectors. 

In conclusion, the perturbations that preserve the PPT property, as well as the ranks (m, n) of p 
and p p , to first order in e, are the solutions of eq. (135T ). 

We may want to perturb in different ways, for example such that Img p' = Img p, but not neces- 
sarily lmg(p') p = Img p p , we only require (p') p and p p to have the same rank. Then the conditions 
on A are that PA = A and QA = 0, or equivalently, 

(I - P + Q)A = . (36) 

4 Product vectors in the kernel 

Assume that p £ V. Recall that the equations pw = and pw = for w 6 are equivalent, and 
so are the equations p p w = and p p w = 0, because p > and p p > 0. Taken together with the 
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identity 

(x ® y)V (u ® v ) = (x <g> v*)^p p (u ® y*) (37) 

this puts strong restrictions on p when we know a number of product vectors in Ker p. 

Assume from now on that w is a product vector, w = u (g> v. Defining w = u (g> v * we have the 
general relation 

pw = w' p p w . (38) 

Assume furthermore that it; € Kerp, this is equivalent to the condition that w G Kerp p . For any 
z € C w we have the condition on p that pw = 0. In particular, when z is an arbitrary product 
vector, z = x ® y, we have the two conditions on p that 

(a; ® yjtp C" ® *>) = (39) 

and 

(x ® v)tp (u ® y) = (x (8) y*)^p p (u ® v*) = . (40) 



Assume that iOj = Uj eg t>j € Ker p for z = 
i, j, we have the following constraints on p, 



1,2, ... ,n. Then for arbitrary values of the indices 

(41) 

(42) 

Ckiij = i (-AfcZij — (Akiij)^) , (43) 



(it* <8> VjY p (u k ® v k ) = (ui <g> v k y p {u k <S> Vj) = . 
Let us introduce matrices 

A-kUj = (u k <8> vi)(ui <S) Vj) ] , 

and Hermitean matrices 



then the constraints on p are of the form 

TripBkkij) = Tr(pCkkij) = Tv(pB kjik ) = Tr(pC kjik ) = . (44) 

Each equation Tv(pB) = with B ^ or Tr(pC) = with C ^ is one real valued constraint. Of 
course, the constraints in eq. (1441 are not all independent, we have for example that Ckkkk = 0- 

5 Rank (4, 4) PPT states in 3 x 3 dimensions 
5.1 The UPB construction of entangled PPT states 

We will review briefly the construction of a rank (4, 4) entangled and extremal PPT state p in 3 x 3 
dimensions from an unextendible orthonormal product basis (a UPB) of Ker p [ 13 ]. The UPB consists 
of five orthonormal product vectors Wi = NiUi <g> vi with the property that there exists no product 
vector orthogonal to all of them. We include real normalization factors iVj here because we want to 
normalize such that wjwj = 5ij without necessarily normalizing the vectors u, and V{. 

The orthogonality of the product vectors Wi follows from the orthogonality relations u\ _L 112 -L 
U3 _L 114 _L U5 _L u\ and v\ _L U3 _L v$ _L V2 -L V4 -L v\. There is the further condition that any three 
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vectors U{ and any three u$ are linearly independent. The five dimensional subspace spanned by these 
product vectors is the kernel of the density matrix 




(45) 



i=l / 

which is proportional to a projection operator. The partial transpose of p is 



with {5j = iVjiij ® t>*. Thus we have both p > and p p > by construction. Note that if all the 
vectors V{ are real, v* = Vi, then p is symmetric under partial transposition, p p = p. 

By a unitary product transformation as in eq. <130T > we may transform the above orthogonal UPB 
into the standard unnormalized form ||9] 

1 d c\ 

1 1 c , (47) 
-c 1 d/ 

with a, b, c, d as positive real parameters. The following quantities determine these parameters, 

det(uiU2ii4) detluiusu^) 2 
s i = — m \ i 77 \ = a , 




det(uiii2U5 


) det(iti«3«4) 


det(uiU2U3 


) det(u2U4Us) 


det(uiU2U4 


) det(u 2 u 3 u 5 ) 


det(viv 2 v^) 


det(viV4Vs) 




det(v 1V3V4) 


det(viV3V5) 


det(v 2 V3V4) 



a 2 



s 2 = —JTZr. ~ ~ \ ~ ~ \ = — > (48) 

and 



s 3 = TT 7 x , 77 7 = C 2 



det(vit)2f3) det(v3U4f5) c 2 



s 4 = ^-TTT! ~ ~ \ TZTF7. 7 7. \ = ■ ( 49 ^ 



These ratios of determinants are invariant under general SL ® SL transformations, as well as indepen- 
dent of the normalization of the vectors. 

In Ref. we presented numerical evidence that every entangled rank (4, 4) PPT state is SL (g) SL 
equivalent to some state of the form of eq. (|45T ) with real product vectors as given in eq. (|47"T ). 

This means that the surface of all rank (4, 4) entangled PPT states has dimension 36. We count 32 
degrees of freedom due to the SL(3, C) <g> SL(3, C) transformations, plus the 4 real SL SL invariant 
parameters a, 6, c, d in eq. (l47l) . 



5.2 A different point of view 

We will present here the construction of an entangled PPT state p of rank (4, 4) as seen from a different 
point of view. When p has rank 4 it means that Kerp has dimension 5. A generic 5 dimensional 
subspace in C 9 = C 3 C 3 has a basis of product vectors. In fact, it contains exactly 6 product 
vectors, any 5 of which are linearly independent. By eq. (1221 . 5 is the limiting dimension for which 
the number of product vectors is nonzero and finite, and the number 6 is consistent with eq. (l23l . 
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Any set of 5 product vectors Wi = m <8> v%, orthogonal or not, may be transformed by an SL ® SL 
transformation to the standard unnormalized form 




10 1 1 
10 1/ 
11s 



(50) 



with p, q, r, s as real or complex parameters. We impose here again the condition that any three m and 
any three Vi should be linearly independent. There is always a 6th product vector which is a linear 
combination of the above 5, 



ps 

1 



q 

r 



qr 

_s 

- s 

p ) 



( 



^6 
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(51) 



The values of the above invariants as functions of the new parameters are 

p r — s r 

si = — , s 2 = q - 1 , s 3 = , s 4 = . (52) 

q s 1 — r 

The parameters p, q, r, s are actually new invariants, they can not be changed by SL ® SL transforma- 
tions. 

If the values of p, q, r, s are such that the invariants si, S2, S3, S4 are all real and strictly positive, 
then we may use eq. (|48T ) and eq. d49l ) to find corresponding values of o, b, c, d, and we may transform 
from the non-orthogonal standard form in eq. (150^ to the orthogonal standard form in eq. (|47V which 
in turn defines the rank (4, 4) state in eq. (l45l) . 

We see from eq. (l52l that the invariants are all strictly positive if and only if p, q, r, s are all real, 
and p<0, q > I, < r < 1, 0<s<r. These inequalities define the regions marked 1 in the (p, q) 
and (r, s) planes as plotted in Fig. [T] 

As discussed in ref. ||9l there are 10 permutations of the product vectors wi for i = 1, 2, . . . , 5 
which preserve the positivity of the invariants. These permutations form a group G which is the 
symmetry group of a regular pentagon, exemplified by the rotation, or cyclic permutation, Wi 1— > Wi 
with 

Wi = W 5 , W 2 = Wi , W3 = W 2 , W4 = U>3 , U>5 = ^4 , (53) 

and the reflection 

wi = w 4 , w 2 = w 3 , w 3 = u> 2 , W4 = wi , w 5 = w 5 . (54) 

For short, we write the rotation as 51234 and the reflection as 43215. 

There are altogether 5! = 120 permutations of the 5 product vectors wi, and they fall into 12 
classes (left cosets of the group G as a subgroup of the permutation group S§) which are not trans- 
formed into each other by G. We number the classes from 1 to 12, and pick one representative from 
each class as follows, 

1 : 12345 2 : 13245 3 : 21345 4 : 23145 5 : 31245 6 : 32145 
7: 12435 8: 14235 9: 21435 10: 24135 11: 13425 12: 14325 
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Figure 1: Regions for the parameters p, q, r, s defined in eq. (l50l such that the product vectors u, <g> V{ 
for i = 1, 2, . . . , 5 are in the kernel of a rank (4, 4) extremal PPT state. The (p, q) plane is divided 
into 12 regions, with 12 corresponding regions in the (r, s) plane. The numbers 1 to 12 refer to the 
permutations of product vectors given in eq. (l55l l. For example, if (p, q) is in region 7, (r, s) must also 
be in region 7. 

Each of these 12 classes defines a positivity region in each of the two parameter planes, where all 4 
invariants si, S2, S3, S4 computed from the permuted product vectors are positive. The 12 regions are 
disjoint and fill the planes completely, as shown in Fig. [TJ On the border lines between the regions the 
condition of linear independence between any three u vectors and any three v vectors is violated. 

To summarize, we have learned how to test whether a set of 5 product vectors wi = Ui<S>Vi, which 
is generic in the sense that any three u vectors are linearly independent and any three v vectors are also 
linearly independent, span the kernel of a rank (4, 4) PPT state. We transform to the standard form 
defined in eq. (1501 ). by a product transformation and normalization. Then the necessary and sufficient 
condition is that the parameters p, q, r, s are all real, and that the parameter pairs (p, q) and (r, s) lie 
in corresponding regions in the parameter planes, as shown in Figure [T] 

It should be stressed that these conclusions are based in part on numerical evidence, and we have 
no analytical proof which is complete in every detail. 

We will discuss next how to reconstruct a PPT state from the product vectors in its kernel. 

5.3 Matrix representation relative to a non-orthonormal basis 

Consider a basis {e-i} consisting of vectors that need not be orthonormal. The scalar products 

9ij = 4 e i ( 56 ) 

define the metric tensor g as a Hermitean matrix. In the usual way, we write the inverse matrix g _1 
with upper indices, so that 

5>%* = 4- ( 57 ) 

3 
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We define the dual vectors e- 7 such that 



(ey = Y,9 l] e}. (58) 

3 3 

They satisfy the orthogonality relations 

( e <)t e . = e t e i = ^ ? (59) 

and the completeness relation 

1 = £ e^ej = £ e>'ej = £ e^e')* . (60) 

Using the dual basis vectors and the completeness relation we may write any matrix A as 

A = ^2e i A ij (e j ) t with Aij = e}Aej. (61) 

i,3 



5.4 Conditions on p from product vectors in Ker p 

It is possible to construct the rank (4, 4) PPT state p directly from 5 product vectors in Ker p without 
transforming first the product vectors to the orthogonal form. We now describe this construction. 

Given three product vectors wi = ui ® Vi in Ker p, with the restriction that all three Ui and all 
three Vi are linearly independent. Then we have the following product basis of C 9 , not necessarily 
orthonormal, 



e ij = u i ®v j , ij = 11,12,13,21,22,23,31,32,33 . (62) 
With respect to this basis we may define matrix elements of p like in eq. (|6TT >. 

PijM = e-ijPe-ki ■ (63) 

In order to count the independent constraints, it is convenient use the standard form of the product 
vectors defined in eq. (1501 . Now all the constraints from eq. (I44l imply that 
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with real diagonal elements oj, ci2, • • • , a>6 and complex off-diagonal elements b%, 62, ■ ■ ■ , b§. This 
Hermitean 9x9 matrix contains 18 real parameters, which means that there are altogether 81 — 18 = 63 
independent real constraints. 
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Including the fourth product vector from eq. (1501 gives additional constraints pw^ = 0, or explic- 
itly written out, 

at + h + b 2 = , 
b\ + a 2 + b 3 = , 
a 3 + 6 4 + 5 5 = , 

63 + 64 + a 4 = , (65) 
b* 5 + a 5 + b e = , 
b 2 + be + a 6 = . 

Here we have simplified slightly by complex conjugating the 4th and the 6th equation. These are com- 
plex equations, to be split into real and imaginary parts. The real parts are 6 independent equations, 
whereas the complex parts are only 5 independent equations. However, we get another independent 
equation as the imaginary part of, for example, the complex equation 

(ui <g) V4)" 1 " p (ui <g> v 2 ) = ai + b\ + b 2 = . (66) 

The end result is that all the off-diagonal matrix elements hi have to be real. Altogether, we get 12 
independent real constraints, 6 from the real parts and 6 from the imaginary parts of the equations. 

Thus, including the 4th product vector in Ker p increases the number of independent real con- 
straints from 63 to 75, and reduces the number of real parameters in p from 18 to 6. 

The generic case with 5 product vectors is that there are 8 1 independent constraints, leaving only 
the trivial solution p = 0. In order to end up with one possible solution for p we have to choose the 
parameters p, q, r, s to be real. 

When we choose real values for p, q, r, s, there is always (generically) exactly one solution for p, 
that is, there are 80 independent constraints. The problem is that this uniquely determined matrix p, 
or its partial transpose, has in general both positive and negative eigenvalues. 

The condition to ensure that both p > and p p > (with the proper choice of sign for p), when 
the parameters p, q, r, s are real, is that the pair (p, q) and the pair (r, s) must lie in corresponding 
parameter regions, as shown in Fig. Q] 

5.5 Separable states of rank (4, 4) 

A separable state of rank 4 has the form 

1 



=1 



with Aj > 0, X)i=i Ai = 1> ipitpi = 1> and ipi = Ci<pi® Xi with C« as a normalization constant. In 
the generic case when any three vectors (pi and any three \i WQ linearly independent, we may perform 
an SL (g) SL transformation and obtain the standard form 




X = I I I I (68! 
In this standard form the x vectors are real, and hence p 1 
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The kernel Kerp consists of the vectors that are orthogonal to all 4 product vectors ipi, and it 
contains exactly 6 product vectors Wi = NiUi® vi, as follows, 

/0 1 1 1 \ / 1 1100 0\ 

u= 1 -1 -1, v = -1 1 10 . (69) 
\ 1 10-10/ V -10-101/ 

Note that the product vectors in the kernel of a separable rank (4, 4) PPT state are not generic, in 
that there are subsets of three linearly dependent vectors both among the u vectors and among the v 
vectors. 

The surface of separable states of rank 4 has dimension 35 = 32 + 3, where 32 is the number of 
parameters of the group SL(3, C) <g> SL(3, C) and 3 is the number of independent coefficients Aj in 
eq. (|67]>. 

6 Rank (5,5) PPT states in 3 x 3 dimensions 
6.1 The surface of rank (5, 5) PPT states 

Since we believe that we understand completely the rank (4, 4) entangled states in dimension 3 x 3, a 
natural next step is to try to understand the (5, 5) states in the same dimension. 

As discussed in the previous section, a generic 5 dimensional subspace in 3 x 3 dimensions con- 
tains exactly 6 product vectors, which can be transformed by SL ® SL transformations, as in eqs. (I30l l 
and (l3Tb . into the standard form given in eqs. (l50l i and (|5TV with SL®SL invariant complex parameters 
p, q, r, s. Thus, each such subspace belongs to an equivalence class under SL <g> SL transformations, 
and the equivalence classes are parametrized by 8 real parameters. There is a discrete ambiguity in 
the parametrization, since it depends on the ordering of the 6 product vectors. 

In one given generic 5 dimensional subspace we may construct a 5 dimensional set of rank (5, 5) 
separable states as convex combinations of the 6 product vectors in the subspace. However, we find 
numerically that the dimension of the surface of rank (5, 5) PPT states with the given subspace as 
image is not 5 but 8. We compute this dimension in the following way. 

We search numerically for one rank (5, 5) state p, for example by the methods described in JSJ. 
The state we find will typically be entangled and extremal in V. From this state p we compute the 
projections P, Q, P, Q as described in Sections |2] and |3j and we look for perturbations p' = p + eA, 
with Tr A = 0, where A satifies both equations PA = PAP = A and QA = (QA P Q) P = 0, or 
equivalently eq. (I36l ). 

(I - P + Q)A = . (70) 

The number of linearly independent solutions for A is the dimension of the surface of rank (5, 5) PPT 
states at the point p. 

We believe that the dimension 8 can be understood as follows. We may fix both subspaces Img p 
and Img p p , this means that we fix the projections P and P and determine p as a solution of the 
equation 

(2I-P-P)p = 0, (71) 

with Tr p = 1. Then there is typically no solution at all for p, solutions exist only for special pairs of 
subspaces. If now the two 5 dimensional subspaces are chosen in such a way that a solution exists, 
then the solution is (typically) unique, and the uniqueness means that p is an extremal point of V. 
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We may fix instead Img p but not Img p p , only the rank of p p . Then there is a set of solutions for 
p described by 8 real parameters. It is a natural guess that the role of these 8 parameters is to specify 
the SL (£> SL equivalence class to which the 5 dimensional subspace Img p p belongs. 

In fact, when we fix Img p there is no degree of freedom left corresponding to SL ® SL transfor- 
mations. This is so because the set of product vectors in Img p is discrete and can not be transformed 
continuously within the fixed subspace Img p. Hence, the only way to vary the subspace Img p p 
without varying Img p is to vary the equivalence class of Img p p . 

Figure [2] shows a two dimensional section through the set of density matrices. The section is 
defined by the maximally mixed state, by a randomly selected rank (5, 5) entangled and extremal PPT 
state p, and by a direction A through p such that the perturbed state p' = p + eA is a rank (5, 5) PPT 
state for infinitesimal positive and negative e, and has Img p' = Img p even for finite e. The figure 
illustrates the fact that the difference between the sets V and S is small. It also illustrates that the 
difference between V and S is largest close to extremal entangled PPT states. 



0.4 




-0.3 1 1 1 1 1 1 1 1 

-0.3 -0.2 -0.1 0.1 0.2 0.3 0.4 

Figure 2: Two dimensional section through T>, the set of density matrices in 3 x 3 dimensions. The boundaries 
of <S, the set of separable states, and of V, the set of PPT states, are both drawn as thick lines. The boundaries 
of T> and of 1> p , the set of partially transposed density matrices, cross at two points, they are drawn as thin lines 
where they do not coincide with the boundary of V. The origin, marked by a small circle, is the maximally 
mixed state. The point marked by a small square is an extremal rank (5, 5) PPT state. The boundary of V is 
drawn through this point as a thin straight line. The boundary of V p is drawn thick at this point, beacause it is 
also the boundary of V . The lightly shaded region around the origin is S. The small and more darkly shaded 
region close to the (5, 5) state is the difference between V and S. Away from this small region, the boundaries 
of V and S are indistinguishable in the plot. 
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If we want to allow both Img p and Img p p to vary, but still require the ranks of p and p p to be 5, 
the equation to be solved for the perturbation A is 



{Q + Q)A = 0. 



(72) 



In this case the number of linearly independent solutions for A is 48, and this is the dimension of the 
surface of all rank (5,5) PPT states. 

We understand the dimension 48 as follows. There are 8 + 8 = 16 parameters for the SL <S> 
SL equivalence classes of the subspaces Img p and Img p p . And there are 32 parameters for the 
SL(3, C) (g> SL(3, C) transformations. 

To summarize, an extremal (and hence entangled) rank (5, 5) PPT state p is uniquely determined 
by eq. (1711 1. as soon as we specify the 5 dimensional subspaces Img p and Img p p . Each subspace 
Img p and Img p p is determined by 8 real SL®SL invariant parameters and an SL(g>SL transformation. 

According to our understanding, based on numerical studies, which is so far only a plausible 
hypothesis, the 8 invariant parameters can be chosen independently for Img p and Imgp p , but the 
SL <g) SL transformations can not be chosen independently. However, we do not know the relation that 
clearly exists between Img p and Img p p . 

In other words, we do not know any explicit procedure for constructing the most general rank 
(5, 5) PPT states. Therefore we turn next to a more restricted problem. 

6.2 Perturbing from rank (4, 4) to rank (5, 5) 

We will see now how to construct (5, 5) states that are innnitesimally close to (4, 4) states. 

Consider once more an infinitesimal perturbation p' = p + eA, this time with p as the rank (4, 4) 
state defined in eq. (143T ). involving the standard real product vectors defined in eq. (l47l . The most 
general case is equivalent to this special case by some SL (g) SL transformation. 

An extra bonus of this special choice of p is that p p = p. In the notation we have used above, we 
have projections P = P on Img p = Img p p and Q = Q on Ker p = Ker p p . 

Conditions on the perturbation matrix A 

By eq. (fl~2l) . the condition for p' to have rank 5 is that 



QAQ = aww 



(73) 



where a is a real number, a ^ 0, and 



5 




(74) 



with complex coefficients Cj such that 



5 




(75) 



Similarly, the condition for (p') p to have rank 5 is that 

QA P Q = QA P Q = f3zJ , 



(76) 
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where (3 is real, (3 ^ 0, and 



^^diWi with z 1 z = |dj| 2 = 1 . 



(77) 



8=1 



1=1 



Note that the possibilities that p' has rank either (4, 5) or (5, 4) are included if we allow either aor/3 
to be zero. 

By eq. (T23T ). there is one extra product vector in Ker p = Ker p p , it may be written as 



wq = aiWi 
i=i 



(78) 



this time with real coefficients dj. Since Wi = NiUi <S> V{ with iVj real and Vi real for i = 1, 2, . . . , 6, 
we have for any Hermitean matrix A that 



wjAwi = w}A F Wi . 



(79) 



By the definition of the projection Q we have that Qwi = Wi for i = 1, 2, . . . , 6. It follows then from 
eq. (|73]» that 



wjAwi = w^QAQwi = a \w}w\ 2 , 

and from eq. (1761 that 

w}A F Wi = w}QA p Qwi = /3 \w}z\ 2 
Together with eq. (1791 this gives the equations 

a \ci\ 2 = (3 \di\ 2 
for % = 1,2,... =5, and the 6th equation 

5 



a 



i=i 



i=i 



(80) 



(81) 



(82) 



(83) 



It follows further that 



J2a\a\ 2 = J2P\di\ 2 = P, 



a = a lc 



i=i 



and that 



|cj| = \di\ for % = 1, 2, . . . , 5 . 



(84) 



(85) 



Thus, the coefficient di can differ from Cj only by a phase factor. The total of 5 phase factors are 
reduced to 4 independent phase factors by the extra equation 



£ 

2 = 1 



(Xi Co 



y^^jdj 



(86) 



i=l 
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For infinitesimal values of e, both p' and (p') p will have four eigenvalues innnitesimally close to 
1/4 and one eigenvalue close to zero, which is ea for p' and e/3 for (p') p . This eigenvalue is the same 
for p' and (p') ', since a = (3. With a > this means that both p' > and (//) P > for e > 0, but 
not for e < 0. Thus, we get automatically a PPT state of rank (5, 5), we never get rank (5, 4) or (4, 5). 
Also it never happens that p' is not a PPT state for the reason that one of p' or (p') p has a negative 
eigenvalue. 

For a more general rank (4, 4) state p, which is obtained by some SL <g) SL transformation from a 
state of the special type discussed here, the smallest positive eigenvalues of p' and (p') p are no longer 
equal. But they are still tied together in such a way that they go to zero simultaneously when we move 
along the surface of (5, 5) states and approach its boundary. The boundary must therefore consist of 
(4, 4) states. 

Computing A 

Define W = ww^ and Z = zz', in the same notation as above. These are both projections, W 2 = W 
and Z 2 = Z, with QW = WQ = W and QZ = ZQ = Z. It follows from eq. <U3 and eq. (5) that 

WAW = WQAQW = aW 3 = aW = QAQ , 

ZA P Z = ZQA P QZ = (3Z 3 = (3Z = QA P Q . (87) 
Like in eq. (fTBl and eq. (1321 we define 

PI = PXP , QX = QXQ , PX = (PX P P) P , QX = (QX P Q) P , (88) 
and furthermore, 

WX = WXW , ZX = (ZX P Z) P . (89) 

We may also define S = Q — W and S = Q — Z, these are again orthogonal projections on H^. 

The least restrictive conditions we may impose on A are now that both S^4 = and S^4 = 0, or 
equivalently, 

(S + S)A = . (90) 

To compute A from this equation we introduce an orthonormal basis in the real Hilbert space H]y. 
Relative to this basis, the operator S + S is represented by a real symmetric positive semidefinite 
matrix, which has a complete set of real eigenvectors with real eigenvalues. We choose A as an 
eigenvector of S + S with eigenvalue zero. 

Apart from the trivial solution A = p, we find 37 linearly independent solutions of eq. (l90l . 36 
out of these 37 are perturbations that give p' = p + eA as a rank (4, 4) state both for e > and e < 0. 
They do not depend on either vector w or z, since they satisfy the conditions Q^4 = and Q^4 = 0. 
But because W = WQ and Z = ZQ they also satisfy the conditions 

WA = WQA = , ZA = ZQA = , (91) 

and hence eq. d90b - The number 36 is the dimension of the surface of rank (4, 4) extremal PPT states, 
as noted in Subsection 15.11 The 37th independent solution is the one giving a rank (5, 5) extremal 
PPT state. 
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A more restricted class of perturbations consists of those where we fix the 5 dimensional subspace 
Img p' to be the direct sum of the 4 dimensional subspace Img p and the one dimensional subspace of 
the vector w. The projection on Img p' is then 



P 5 = P + W , 



(92) 



and the partial condition on A is that P$A 



A, when we define 



P 5 X 



(93) 



The full condition on A is that 



(P 5 



S)A = A. 



(94) 



Again apart from the trivial solution A = p, we find 5 linearly independent solutions of eq. (l94l . 
of which 4 give p' as a rank (4, 4) state both for e > and e < 0. The 5th independent solution is the 
one giving p' as a rank (5, 5) extremal PPT state. 

The 4 directions that only give new (4, 4) states are easily identified, since they do not depend on 
the vector z. To find them we simply repeat the calculation with a "wrong" z, violating the condi- 
tions (|85T ) and (l86l ). In this way we find no (5, 5) state, but we find the same set of perturbations into 
(4, 4) states. The number 4 is the dimension of the surface of (4, 4) states with image within the fixed 
5 dimensional subspace projected out by the projector P5. 

There is a natural explanation of why this surface has dimension 4. In fact, when we fix P5 and 
look for (4, 4) states with image within this fixed 5 dimensional subspace, we eliminate all degrees of 
freedom corresponding to SL (g> SL transformations. But we still allow variations of the 4 real SL ® SL 
invariant parameters that are needed to define a rank (4, 4) state. 

We conclude that for fixed vectors w and z there is one direction away from the surface of rank 
(4, 4) extremal PPT states and into the surface of rank (5, 5) extremal PPT states. 

For a fixed vector w there is a 4 parameter family of acceptable vectors z. Recall that these 4 pa- 
rameters determine the 5 relative phases between the coefficients a in eq. (1741 and the corresponding 
coefficients di in eq. (I77T ). 

The vector w is an arbitrary vector in the 5 dimensional kernel of the unperturbed state p, hence 
it contains 4 complex parameters, or 8 real parameters, after we take out an uninteresting complex 
normalization factor. Altogether, there are 8 + 4 = 12 independent directions away from the 36 
dimensional surface of rank (4, 4) PPT states and into the surface of rank (5, 5) PPT states. 

When we perturb an arbitrary rank (5, 5) PPT state in such a way that we preserve the ranks of 
the state and its partial transpose, we find numerically that the surface of rank (5, 5) PPT states has 
dimension 48. The fact that 48 = 36 + 12 is consistent with the hypothesis that we can reach every 
rank (5, 5) PPT state if we start from a rank (4, 4) PPT state and move continuously along the surface 
of rank (5,5) PPT states. 

7 Rank (6,6) entangled PPT states in 4 x 4 dimensions 

We will discuss in some detail one more example of the relation between PPT states and product 
vectors. According to eq. d22l . the rank (6, 6) PPT states in 4 x 4 dimensions represent just the 
limiting case with a finite number of product vectors in the kernel, in this respect they are similar to 
the rank (4, 4) states in 3 x 3 dimensions. 
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The kernel of a rank 6 state in 16 dimensions has dimension 10, and the generic case, according 
to eq. (|23T ). is that it contains exactly 20 product vectors, any 10 of which are linearly independent. 
We will see here that the product vectors in the kernel put such strong restrictions on the state that the 
rank (6, 6) PPT state may be reconstructed uniquely from only 7 product vectors in its kernel. 

To see how it works, take a set of product vectors in 4 x 4 dimensions. We may take random 
product vectors, or else a set of product vectors with the special property that they belong to Ker p 
where p is a rank (6, 6) PPT state. We find numerically that the number of constraints generated by 
fewer than 7 product vectors is the same in both cases. We find the following numbers. 

From 4 product vectors assumed to lie in Ker p for an unknown p, or actually lying in Ker p for 
a known p, we get 172 independent constraints on p of the form given in eq. (I44l . These constraints 
leave 84 free real parameters in p, before we normalize and set TV p = 1. 

From 5 product vectors in Ker p we get 205 independent constraints, leaving 51 parameters in p. 

From 6 product vectors in Ker p we get 234 independent constraints, and 22 parameters in p. 

Finally, 7 product vectors in Ker p give either 255 or 256 independent constraints, and either 1 or 
real parameters in p. If there is one parameter left, it is a proportionality constant, to be fixed by the 
normalization condition Trp = 1. 

The standard form of 7 product vectors in 4 x 4 dimensions, generalizing eq. (I50l . is the following, 



/I 











1 


1 


M 




/I 











1 


1 


1 \ 





1 








1 


Pi 


Pa 







1 








1 


P7 


Pio 








1 





1 


P2 


P5 


v = 








1 





1 


PS 


Pn 


\o 








1 


1 


P3 


PeJ 




Vo 








1 


1 


P9 


P12/ 



(95) 



There are 12 complex parameters pi,p2, ■ ■ ■ ,Pu, that is, 24 real parameters. These are invariant in 
the sense that we can not change them by SL(4, C) <g> SL(4, C) transformations. 

Not just any arbitrary set of 7 product vectors defines a rank (6, 6) PPT state. We arrive at this 
conclusion not only because we find numerically that 7 generic product vectors allow only p = as 
solution of all the constraint equations, but also because a dimension counting shows that we need 
less than 24 invariant parameters in order to parametrize the rank (6, 6) PPT states. 

Take one known rank (6,6) PPT state p and perturb it into another rank (6,6) PPT state p' = 
p + eA with e infinitesimal. Here A must be a solution of eq. d35l) . with operators Q and Q defined 
relative to p as explained. The number of linearly independent solutions for A, found numerically, is 
76, including the trivial solution A = p. This shows that the surface of rank (6, 6) PPT states has 75 
real dimensions. 

Of these 75 dimensions, 60 dimensions result from product transformations p h-> VpV^ with 
V = Va <8> Vb and Va, Vb € SL(4, C). The remaining 15 dimensions must correspond to 15 SL ® SL 
invariant parameters of the 7 product vectors. 

It is also worth noting that, by the counting explained in the next section, the set of 6 dimensional 
subspaces of C 16 has real dimension 16 2 — 6 2 — 10 2 = 120, much larger than the dimension 75 of 
the surface of (6, 6) states. Thus, not every 6 dimensional subspace of (D 16 is the host of a rank (6, 6) 
PPT state, as one would expect from the analogy to the case of the rank (4, 4) PPT states in C 9 . 

8 Dimension counting 

We will describe in this section how to compute numerically the dimensions of surfaces of PPT states 
of given ranks. We list some numerical results, and discuss how they may be understood in most cases 
by a simple counting of constraints, assuming the constraints to be independent. 
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We start with a useful exercise. We want to compute the real (as opposed to complex) dimension 
of the set of all r dimensional subspaces of an N dimensional complex Hilbert space. 

First note that the unitary group U(fc) has k 2 real dimensions. Take an orthonormal basis of the 
Hilbert space. The first r basis vectors define an r dimensional subspace, the orthogonal complement 
of which is defined by the last s = N — r basis vectors. A U(iV) transformation transforms this basis 
into another orthonormal basis, but the U(r) transformations within the first r basis vectors, and the 
U(s) transfomations within the last s basis vectors, do not change either subspace. It follows that the 
dimension of the set of r dimensional subspaces, equal to the dimension of the set of s dimensional 
subspaces, is 

d = N 2 - r 2 - s 2 = 2rs . (96) 

Assuming that we have found a PPT state p of rank (m, n), it lies on a surface of rank (m, n) PPT 
states. We compute the dimension of the surface at this point by counting the number of independent 
solutions A of eq. (f35l) . 

(Q + Q)A = , (97) 

equivalent to the two equations 

Q.4 = QAQ = , QA = (QA P Q) P = . (98) 

We have to throw away the trivial solution A = p. We get a lower bound for the dimension if we 
assume that the constraints on A from the two equations in eq. (l98l l are independent. The equation 
QAQ = represents (N — m) 2 real constraints, since Q is the orthogonal projection on the N — m 
dimensional subspace Ker p. Similarly, the equation QA P Q = represents (N — n) 2 real constraints, 
since Q is the orthogonal projection on the N — n dimensional subspace Ker/5 P . Because the con- 
straints are not necessarily independent, we get the following lower bound for the dimension, 

d> N 2 — (N — m) 2 — (N — n) 2 — 1 . (99) 

Take A r = 3x3 = 9asan example. We find numerically that eq. (1991 holds with equality for 
all ranks from the full rank (m, n) = (9, 9) down to (m, n) = (5, 5). In particular, for rank (5, 5) the 
dimension of the surface is 

d = 9 2 - 4 2 - 4 2 - 1 = 48 . (100) 

By eq. (l96l) the set of 5 dimensional subspaces has dimension 40, hence we should expect to find an 
8 dimensional surface of rank (5, 5) PPT states in every 5 dimensional subspace. And that is actually 
what we find. 

For rank (4, 4) the constraints are not all independent, and we have the strict inequality 

d = 36 > 9 2 - 5 2 - 5 2 - 1 = 30 . (101) 

The set of 4 dimensional subspaces has again dimension 40, hence there can not exist rank (4, 4) 
PPT states in every 4 dimensional subspace. There are 40 — 36 = 4 constraints restricting the 4 
dimensional subspaces supporting rank (4, 4) PPT states, and in each 4 dimensional subspace there 
can exist at most one unique such state. The 4 constraints are the conditions that the 4 parameters 
a, b, c, d in eq. (|4?1 . or p, q, r, s in eq. (I50l . have to be real. 
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If we want to compute the dimension of the surface of rank (m, n) PPT states with fixed image 
space, we have to count the independent solutions of eq. (I36l ). 

(I - P + Q)A = , (102) 

equivalent to the two equations 

PA = PAP = A , QA = (QA P Q) P = . (103) 

The equation PA = A leaves m 2 real parameters in A and represents N 2 — m 2 real constraints, as is 
visualized in eq. (fT71) . The lower bound on the dimension is therefore 

d> N 2 — (N 2 — m 2 ) - (N - nf - 1 = m 2 - (N - n) 2 - 1 . (104) 

In the above example with N = 9 and (m, n) = (5, 5) we find numerically d = 8, as already 
mentioned, so that the inequality in eq. (11041) holds as an equality. With (m, n) = (4, 4), on the other 
hand, we get 

d = > 4 2 - 5 2 - 1 = -10 . (105) 



9 Numerical integration 

In this section we will describe a numerical method for tracing curves on a surface of PPT states of 
fixed ranks (m, n). This is a tool for studying the geometry of the surface, for example by tracing 
geodesies to see how they curve, or studying how the surface approaches a boundary consisting of 
states of lower ranks. 

9.1 Equations of motion 

The perturbation expansion p(t + e) = p(t) + eA for p = p(t) is equivalent to the differential equation 

p = A. (106) 

We use the notation 

p= Tt' p ="dT- (107) 

We defined the pseudoinverse p + in eq. (01), in order to define P = p + p = pp + and Q = 1 — P, 
the orthogonal projections on Img p and Ker p, respectively. There are similar relations for P, the 
projection on Img p p , and Q = 1 — P, the projection on Ker p p . We defined orthogonal projections 
on Hn, the space of Hermitean matrices, in eq. ([TBI and eq. (1321 . 

If X is a Hermitean matrix with Img X C Img p then X = PX = XP, or equivalently, QX = 
XQ = 0. Assuming that these relations hold at any "time" t we may differentiate and get that 

X = PX + PX = XP + XP . (108) 

Equivalently, 

QX = PX , XQ = XP . (109) 
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Multiplication by Q from the left and from the right gives that 

QXQ = 0. (110) 

It follows further that 

X = (P + Q)X(P + Q) = PXP+XP+PX . (Ill) 

The special case X = p gives the equation 

QAQ = (112) 

as a consistency condition for eq. (11061 ) with the relations p = Pp = pP. Eq. (11121 ) is the same as 
eq. (ITTb . the condition for the rank of p to be constant. We may want to replace it with the stronger 
condition that Img p should be constant, eq. ©, 

PAP = A. (113) 

Setting X = p in eq. d 109b gives the equations 

QA = Pp, AQ = pP, (114) 

and multiplication by p + gives that 

QAp + = PP, p + AQ = PP. (115) 

Differentiating the equation P = P 2 gives that P = PP + PP, hence 

P = QAp + + p + AQ . (116) 

Differentiating the relation p + = p + pp + we get that 

p+ = p+p + p+Ap + + Pp + . (117) 

When we left and right multiply here by P we obtain the relation 

Pp + P = -p + Ap + . (118) 

Hence, using eq. (11111) with X = p + , together with eq. (11161 ). we get that 

p + = QA(p+) 2 + (p + ) 2 AQ - p+Ap + . (119) 

The equations (1106 1 ). (II 16I ). and dll9l ) may be integrated together, as soon as we specify how to 
calculate A as a function of p. There are, of course, equations similar to (11161 and (11191 that hold for 
the projection P related to the partial transpose p p , and for the pseudoinverse (p p ) + . 

As a specific example, consider how to generate a curve p = pit) lying on the 48 dimensional 
surface in of rank (5, 5) PPT states in 3 x 3 dimensions passing through a given state p(0). We 
then have to satisfy the two conditions on A that 

Q.4 = QAQ = , QA = (QA P Q) P = . (120) 

Or equivalently eq. (135T ). 

(Q + Q)^ = 0. (121) 

Alternatively, we may want to generate a curve that follows the 8 dimensional surface in Hn of 
rank (5, 5) PPT states such that the 5 dimensional subspace Img p is kept fixed, but the 5 dimensional 
subspace Img p p is allowed to change. This means that we replace the condition QAQ = by the 
condition PAP = A. The single condition to be satisfied is then eq. (l36l ). 

(I - P + Q)A = . (122) 

In this case P is constant but Q = Q(t) may vary as a function of t. 
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9.2 Geodesic equations 



Both conditions (11211 ) and (11221) are of the form 



TA = . 



(123) 



Differentiating this equation gives that 



TA + TA = . 



(124) 



It follows that 



A = B - T + TA , 



(125) 



where T + is the pseudoinverse of T, and B is an arbitrary Hermitean matrix with TB = 0. 

By definition, a geodesic on an embedded surface (think of a great circle on the surface of a 
sphere as an example) is a curve which does not change its direction on the surface. Hence, it changes 
direction in the embedding space only as much as it has to in order to stay on the surface. This would 
mean that we choose B = in eq. (1125 1 ). Or if we normalize A to unit length, fixing Tr A 2 = 1, we 
set B = a A and choose a such that Tr(AA) = 0. 

9.3 Numerical results 

We have done numerical integrations by a standard fourth order Runge-Kutta method. With p and A 
of order one and time steps of order 1CP 4 this gives a precision of order 1CT 16 , which is the machine 
precision. 

Figure [3] shows a geodesic curve pit) on the 8 dimensional curved surface of rank (5,5) PPT 
states with lmgp(t) constant. Figure |4] shows the 5 nonzero eigenvalues of p and p p . The condition 
that one eigenvalue of either p or p p goes to zero defines the boundary of the surface. We see that the 
curve approaches the boundary twice, but turns around each time and continues in the interior. The 
eigenvalue spectra of p and p p are remarkably similar, yet they are not identical. When both p and p p 
simultaneously get one dominant eigenvalue, we interpret it as an indication that p approaches a pure 
product state. 

It is quite natural that a geodesic chosen at random will not hit the boundary, since the boundary 
consists of rank (4, 4) PPT states and has dimension 4, while the surface itself has dimension 8. In 
order to hit the boundary we can not follow a geodesic, we have to integrate the equation p = A and 
choose the direction A in such a way that the smallest positive eigenvalue of p goes to zero. When we 
do so, the smallest eigenvalue of p p goes to zero simultaneously with the eigenvalue of p, although 
the ratio between the two eigenvalues goes to a value different from one. Hence the curve ends at a 
(4, 4) state on the boundary. The explanation for this coupling of eigenvalues of p and p p was given 
in Subsection 16.21 

10 Summary 

The work presented here is part of an ongoing programme to study quantum entanglement in mixed 
states. We have studied here low rank entangled PPT states using perturbation theory and the close 
relation between PPT states and product vectors. 

One result obtained is an understanding of how to construct rank (5, 5) PPT states in 3 x 3 di- 
mensions by perturbing rank (4, 4) states. We use perturbation theory to study surfaces of PPT states 
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of given ranks, and in particular to compute dimensions of such surfaces, for example the surface of 
(5, 5) states. However, it is still an unsolved problem how to construct general rank (5, 5) PPT states 
that are not close to rank (4, 4) states. We are even farther from a full understanding of higher rank 
PPT states in 3 x 3 dimensions, or in higher dimensions. 

A special class of PPT states are those of special ranks so that their kernel is spanned by product 
vectors and contains a finite number of product vectors. We have shown that these states may be 
reconstructed uniquely from a subset of the product vectors in the kernel, and the number of product 
vectors needed may be smaller than the dimension of the kernel. This result raises new interesting 
questions to be answered by future research, for example, how to identify finite sets of product vectors 
that define PPT states with these product vectors in their kernel. 
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Figure 3: A projection of a geodesic curve on the 8 dimensional surface of rank (5, 5) PPT states with a fixed 
image space. We have made a principal component analysis and plotted the two largest principal components. 
The curve starts middle right and ends lower left. 




0.5 1 1.5 2 2.5 3 



Figure 4: Variation along the curve in Fig. |3]of the 5 nonzero eigenvalues of the density matrix (full drawn 
lines) and its partial transpose (broken lines). The abscissa is the arc length along the curve. 
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